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^^ I New supersymmetric partners of tlie modified Poschl- Teller and the Dirac's delta 

^^ , well potentials are constructed in closed form. The resulting one-parametric potentials 

are shown to be interrelated by a limiting process. The range of values of the param- 
eters for which these potentials are free of singularities is exactly determined. The 
construction of higher order supersymmetric partner potentials is also investigated. 
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1 Introduction 

The modified Poschl- Teller potential V{a,x) = —Uo{cosh ax)~'^ , is one of the few exactly 
solvable potentials in quantum mechanics. It was first analyzed by Rosen and Morse Q, 
who found its energy eigenvalues and eigenfunctions. (This potential V{a, x) is indeed, 
an hyperbolic version of what is commonly known as the Poschl- Teller potential [|[.) The 
system is characterized by a finite number of bound states whose spectrum depends on 
the parameters C/o,a > 0, plus a continuun of scattering states (see in |y] an elegant 
group approach); the normalization of the wave functions has been obtained by Nieto Q. 
Barut et al Q studied a three dimensional version of the problem (with an additional term 
proportional to (sinhax)"'^ in the potential), and using the Infeld-Hull factorization Q 
they constructed ladder operators to determine the bound and scattering states from the 
matrix elements of group representations. It is well known that the modified Poschl- Teller 
potential appears in the solitary wave solutions of the Korteweg-de Vries equation Q, and 
also that it can be obtained from supersymmetric (susy) quantum mechanics [^ as the 
susy partner of the free particle potential [^, being a non trivial example of unbroken 
supersymmetry |1C]. 



This potential can be considered as a type of short range potentials, almost vanishing 
over most of their domain, except near zero (where the source resides). The extreme 
case of such short range potentials is the Dirac delta well V{x) = —g5{x), g > 0. Delta 
potentials have been long used in field theories, where the main problem arises from the 
regularization and renormalization of the values for the physical observables predicted by 
the theory itself |11]. The effects of adding a delta function potential on the states of 



a previously known potential have been computed exactly by Atkinson |12| for several 



physical systems. The one dimensional problem of N particles interacting by means of 



delta potentials is one of the simplest many-body problems solved exactly (see |13| and 



references quoted therein). As the delta potential 5{x) is a distribution rather than a 



function, it can be approximated by different families of fmictions, one of them {V{a, x) 
with an appropriate choice of the depth Uq) will be investigated in this paper. 

Concerning the study of new exactly solvable problems in quantum mechanics, in the 
last years there has been remarkable progress along different lines: Darboux transforma- 
tion jl^, Infeld-Hull factorization Q, Mielnik factorization |[T5| , 16, |l^, susy quantum 



mechanics |8|, and inverse scattering theory |18|| , among others. It is worth stressing that 
all of them can be embraced in an elegant algebraic approach named intertwining tech- 
nique p^], which has been successfully applied in the construction of higher order susy 



partners |2C, 21]. The usefulness of the intertwining has been also proved in the study and 
interpretation of black-hole perturbations in general relativity [£^] . 

The main purpose of this paper is to analyze the susy partners associated with the 
modified Poschl- Teller and Dirac delta potentials by using the intertwining technique. In 
Section 2, after a short review of results concerning both potentials, the properties of the 
latter are straighforwardly obtained as a limiting case of the former. In Section 3 we will 
determine closed expressions for two new families of susy partner potentials of Poschl- 
Teller, paying attention to the appearence of singularities (in that case the isospectral 
properties should be considered under the point of view of ||2^). Besides, the limit already 
mentioned in Section 2 is carefully studied here for these two families, obtaining very 
different behaviour for each of them. A direct analysis of the susy delta potential is 
also carried out to check the validity of the previous limits. Finally, in Section 4 other 
varieties of Poschl- Teller susy potentials are shown to be easily obtained using higher order 
intertwining. 



2 The essentials of modified Poschl- Teller potential 

Let us consider the well known one-dimensional two-parametric modified Poschl- Teller 
potential [^], written in the following equivalent forms: 

.r, X ^0 ^^ 2 A(A-l) ga 

V[a,x) = 2 = ~T, — ^ 2 ~ 2 ' « > U, (1) 

cosh ax 2m cosh ax 2 cosh ax 

Henceforth, in all of the three versions, the following conditions are imposed in order 
to have an attractive potential: ^/q > 0, A > 1 or g > 0. We also take for simplicity 
fi /2m = 1, and from the last equality the parameters a, A and g are related by 



The bound states {E < 0) for this potential can be obtained by using the traditional recipe 
of transforming the stationary Schrodinger equation into an hyper geometric equation, with 
parameters 

The general solution is 

V'(x) = (cosh ax) A 2-F1 (a, 6; 1/2; — sinh^ ax) 

+B (sinh ax) 2^1 (a -h 1/2, 6 -M/2; 3/2; - sinh^ ax)\ . (4) 

The normalization condition on these eigenfunctions allows to determine the energy spec- 
trum, which is found to be 

^„ = -a2(A-l-n)^ nGN, 0<n<A-l, (5) 

or just in terms of a and g: 



From (|5|) and A > 1, the energy for n = always belongs to the spectrum of V{a,x). 
Calling A^ the biggest possible value of n in (|6|), the total number of bound states is A^-|- 1. 



From there, the following relationship between the parameters a and A^ holds 

< a < 



(2A^ + 3)2-1 - (2A^ + 1)2-1' 

If A^ = 0, then a > g/A, and there is just one bound state such that it has the lowest 
energy Eq, and for which we have 

ao = co = l/2, 60 = A- 1/2, A / 0, 5 = 0, ^0 = -a^ (A - 1)^ (7) 

The corresponding normalized wave function turns out to be 



^A_../•l^ 1.1. ...2 \_ «r(A-l/2),_^_,,„, 



■00(2;) = Co (cosh ax) 2i^i( -,A- -;-;-sinh axj=J (cosh ax) . (8) 

It is also interesting to remark that (Q) is a transparent potential (i.e. the reflexion 
coefficient is equal to zero) when the following condition is verified |2^: 



5 = 2aA;(A; + l), A; = 0,1,2,... (9) 

Using our notation this implies 

X = k + 1, k = l,2,... (10) 

When A; = we have the case of the free particle, which was excluded from the very 
begining. 

The last form of the potential V{a, x) given in (|l]) can be directly related to the Dirac 
delta potential Voi^) = —g5{x) just by taking the limit a — > 00: 

lim Via^x) = —g5{x). (11) 

a — ^00 

As it is well known, this J-well has a unique bound state with eigenvalue 

Es = -{g/2)\ (12) 

and normalized wave function 

V'5(x) = v^ e-3\^'\l\ (13) 



These results can be obtained directly from the modified Poschl- Teller potential by proving 
that the limiting relationship (|ll| ) between both types of potentials is also inherited by 
their eigenfunctions and energy eigenvalues. Indeed, from the the Poschl- Teller ground 
state energy level (|7|) we have 

limi^o = -lim^f^-4 + ---V = -f-V=^^- (14) 

A similar analysis can be done for the ground state eigenfunction (|8|). Observe that 
■00(0) = Co; if X 7^ 0, then the limit a — > cxd is equivalent to ax — > ±00, hence it can be 
calculated from the asymptotic form of the function. To find it, we use 

g»\x\ ^a\x\ 

coshax ~ , sinhax ~ ± , (a > 0), (15) 



and also the fact that, for 2; ^ —00, we have |2^ 

.F,(„,.;c;.i~r(c){^Ii^(-.r + ^M_(-.)-}, |arg(-.)| < .. (16) 
Then, from (M) and (16), we have for |2:| ^ oo 

,f,{^Im;\;-z')-z~'''^. (17) 

In addition we have the following asymptotic behaviour for a -^ oo: 26o ~ 1 + g/a. Hence 
lim Mx) = Coo lim 2^«-5 g-kl^C^o-i) = c^o e'^l^l/^, (18) 



where Cq -^ Coo • From (^) it is clear that Coo = vW^i and then equation (|lq) states the 
connection between (|^ and ([l3| ) when the limiting relationship (|l^) is satisfied. 

Finally, it is also well known that the 5- well potential is transparent. This property 
is also deduced from the Poschl- Teller potential by taking the limit a ^ co in ([10[), and 
using (pi), which gives the solution A; = 0. 

3 Susy modified Poschl- Teller potential 

Let us consider now the problem of finding the supersymmetric partner of the modified 
Poschl- Teller potential V{a, x). We look for a first order differential operator A = ^+/3(2;) 



and a partner potential V{a, x) such that the fohowing interwining relationship holds: 



■d^ + ^("'"^ 



A = A 



"^ + ^("'^) 



(19) 



The new potential V{a,x) is related to V{a,x) through the following susy relationship 



V{a,x) = V{a,x) +2p'{x), 



(20) 



where f3{x) is a solution of the Riccati equation 



p^{x)-p'{x) = V{a,x)-e, 



(21) 



with e an integration constant, which turns out to be the factorization energy. There is 
an immediate particular solution of equation ( [2l| ) in the form of an hyperbolic tangent, 
Po = D tanhax, with D depending on a. The introduction of /?o in (pi]) gives 



D- 



-aX, D = — a(l — A). 



(22) 



Therefore, we have two different particular solutions of (21) 



/3q (a, x) = D tanh ax, 



(23) 



associated with two different factorization energies 



..._,„*,.. _^,,i.y,,M 



(24) 



Remark that these factorization energies can be formally identified with two values of the 
spectrum formula (^): 

e- = Eo, e+ = E^i. (25) 



Then, the general solutions of the Riccati equation (21) for the above factorization energies 
can be found to be 

/?+ (a, x) = D+ tanh ax-^ln(l-C T (cosh ay)^^^/" dy) , (26) 

f3^{a,x) = D- tanh ax - — In ('l - ^ r {cosh ay f^~ ' '^ dy\ , (27) 



where (, ^ are two new independent integration constants; when they are taken to be zero, 
we recover the particular solutions /3q (x). It must be clear that, in fact, we have obtained 
two different families of intertwining operators 

A^ = ^ + P^{x\ Aj = -^ + P^{x), (28) 

generating two different families of susy partners (pO) of the potential (1^). Up to now, it 
has been usual to consider only the susy partners of a given potential constructed by taking 
particular solutions of the Riccati equation |^, |^. For the potential we are dealing with, 
the cases ( = or S, = give interesting results, and will be obtained just as byproducts 
of (|2q)-(p7|). In principle, it is possible to find solutions associated to other factorization 
energies, but they will produce very awkward expressions, without adding new relevant 
information to the problem we are dealing with. 

Let us analyze now the results coming out when the general solutions (p6|)-(|27|) are 
taken into account. Observe that the integrals appearing there can be expressed in a 
closed form as 

/ (cosh ayydy = ai^i -^, -g; 1 - ^; -e^""^ + constant. (29) 

J aq \ 2 2 J 

This expression will be used next to construct the closed form of the susy partner poten- 
tials. 

3.1 The two parametric family of potentials V^{a,x) 



In this case, taking into account (22) and (|2q), we see that the exponent in ( |29| ) is negative, 
indeed q = — 2A. The definite integral exist in the whole real axis, and in this case we can 
define the function 

2ar(A + l/2)' 
(30) 



M(A, a,x) = / (cosh ay)'"^^ dy = ^ — 2F1 (A, 2A; 1 + A; -e^ 

7o 2aA V 



A typical plot of M(A, a, x) is shown in Figure 1 for general values of the parameters a 
and A. It is quite clear that this function is odd in the variable x, and it is monotonically 



increasing from its minimum value M(A, a, — oo) = — y^r(A)/(2a r(A + 1/2)) to its max- 
imum value M(A, a, +00) = |M(A,a, — oo)|. It is interesting to remark the resemblance 
between the form of this function and that of the error function erf (x) = —j= Jq e~^ dy. 
We will immediately see that, indeed, for the Poschl- Teller problem, M(A, a, x) plays 
a role completely analogous to that played by the error function when determining the 



one-parametric family of susy partners of the harmonic oscillator potential [15|. 

If we introduce now the function 

^4-., N (coshax)"^'*' 1 (i , , ^,,/^ ^ , ^ 

"?<^'"-^-> ■= l-(M(X.a..) = "^'"'^ -CM(X,a..>] (31) 

we have for (3'^ in (26) the following expression: 

f5t(a,x) = —aX ta,nh. ax + C^f(X, a, x) = -—In (cosh ax) f!.f(X,a,x) . (32) 

^ ^ ax L "^ J 

From here, and using (|20|), we can evaluate the associated susy partner potential, which 
turns out to be 



V+(a,x) = -a' l + ^ + Jl + ^^ ^ 



*= ' \ 2a V a J cosh^ ax 

}t{X,a,x) tanhax + 2((^ fit 

It is obvious that the singularities of V^ correspond to the singular points of the function 



-4AaCf^^(A,a,x) tanhax + 2(C f^^(A,a,a;))2. (33) 



Q,f{X,a,x). It can be proved that this function is free of singularities in the following 
range of values of C 

'^' M(A,a,+oo) 0Fr(A) ^ ^ 

As we pointed out before, if we compare our results with the pioneering Mielnik's work on 
the harmonic oscillator [15|, one can appreciate that the roles played there by the error 



function and his parameter 7, are performed here by M(A, a, x) and the inverse of C- The 
characteristic features of M(A,a,x) determine the existence of susy partner potentials 
wich are free of singularities, and are therefore almost isospectral to the modified Poschl- 
Teller potential. (A similar analysis can be done for an equivalent integral appearing in 



Let us work now in this range of values of the parameter (^. The potential (^) corre- 
sponds to the following family of almost isospectral Hamiltonians 

H+ ■■=-^ + ^c^(«, ^) = 4 (4)^ + e+. (35) 

It is well known that the eigenfunctions of H7' can be constructed by acting with the 
operator At of (28) on the eigenfunctions tpn of H, which is factorized as 



H = -^ + Via,x) = {Ap^A+ + e+, (36) 

and are given by i/^niC^ x) o^ Af ipnix), provided that Af ipn{x) ^ and V'n (d ^) S L^(M). 
There is also the possibility of an extra eigenfunction ip^{C,,x) satisfying 

(^+)t^+(C,x)=0, (37) 

which will be called "missing state". 

The first point we want to stress is that according to ( p5[) e+ = -B-i, wich is an energy 



level not allowed in the spectrum of the initial Hamiltonian H. Hence, from (3f:) it is clear 
that there is no eigenfunction ip oi H annihilated by Af. Therefore, the eigenfunctions of 
Ht are given by the normalized functions 

i;+{C,x) = {En-e+)-y^A+iPn{x), n = 0,l,..., (38) 



plus the missing state solving (37), which properly normalized reads 



ip+iC,x) = J \ ^ ' ' ^ coshax^ 17+ A,a,x . 
y 2M[X,a,+oo) ^ 



(39) 



Remark that the non-singularity condition (^) appears here again, although in this case 
it is required for the missing state to be normalizable. Note also from (pEj) that ^"'"(C,x) 
is clearly the eigenfunction of Ht with eigenvalue e+. This is the reason why we named 
it missing state. The spectrum of Ht is given by the set {En, n = 0, 1, . . .} plus a new 
level at e+ = E^i. The conclusion is immediate: the family of potentials Vt~{a,x) is not 
strictly isospectral to its susy partner V{a,x): it has the same levels plus an additional 

10 



one which is placed below all of them. Let us remark that, due to the annihilation of 
the missing state (f^{(,x) in (^) by the intertwiner At, the missing state has no susy 
partner, and the couple of almost isospectral Hamiltonians H and Ht corresponds to a 
case of unbroken susy. 

In Figure 2 we plotted the asymmetric double well corresponding to the susy partner 
potential V^{a,x) given by (33), with a = 0.1, A = 3, and ( = 0.0937. The three bound 



states of this potential are also represented with dotted horizontal lines. Note that the 
potentials V^{a, x) present features that makes them interesting for physical applications: 
(i) we are able to know their spectra in an exact form, and by adjusting the parameters 
we can have the desired number of bound states, and (ii) the shape of the potential can 
be modified to allow interesting tunnelling effects. 



Let us consider now the limit a — > oo of the potentials (|33|). There are two terms 
containing the function i7t(A, a,x); using ([T5|), (16), and (^), it can be easily proved 



that it has the following behaviour for large values of a: 

"?(A,„,.)„~{^^";*'- ^f»; (40) 

giving a discontinuous function in the limit, but which is zero almost everywhere. In 
addition, the product aQ,t{X,a,x) — > 4:6{x), and a^f {X, a, x)tanh. ax —>■ 0. Hence, the 
only relevant part in this potential would be that coming from the first term. But it 
diverges very badly as —4:a5{x), and therefore we do not end with a physically interesting 
potential. 

3.2 The two parametric family of potentials Vt~{a,x) 



The study of the other family of potentials coming out from (27) can be done according 
to the lines already followed in the previous subsection. Nevertheless, there are important 
differences between the results obtained in both cases. First of all, from equations (p^) 



and (pTD, the exponent in (29) is now positive: q = 2(A — 1). As a consequence, if we try 



to evaluate the integral in the whole real axis we get a divergent result. Nevertheless, it 
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is useful to introduce a function similar to M{X,a,x), let us call it 

L{X, a, x) = r {cosh ay f'-^-^^dy (41) 

Jo 

^-2aiX-l). V?r(2-A) 

- -^^(^=^)MA^ ' ' ^ ' ' ' J + 2a(A-l)r(i-A)- 

This function is also odd and takes arbitrary positive values for x > and arbitrary 

negative values for x < 0. Using the asymptotic behaviour of the hyper geometric functions 

(|l6|), it is very easy to prove that the limit a — > cxd of ( pl|) is 

/gSl^l - l\ 
lim L(X,a,x) = sgnx, (42) 

Q^OO \ 9 J 

where sgnx denotes the function sign of x, and we have used the fact that a, A and g are 
related through equation (^). We will use L{X,a,x) to define the function 

^ /, X (coshaj;)^^'*'^^) Id,. ^^.^ ^ ,. . 

»,-(A.a.x),= ;_^^J^^^^,^ =---ln[l-a(A.a.x)], (43) 



from which the following expression for /3c in (27) is obtained 



d , , , ,i_A 



/3c (a,x) = a(A — 1) tanhaa; + .^ ilt (A, a, j;) = — In (cosh ax) ^c{X,a,x) . (44) 
Using this expression and (SO) we compute the new susy partner potential 



^ \ 2a \ a cosh ax 



+4a(A- 1)^ n^{X,a,x) tanhax + 2(^ n'^{X,a,x)). (45) 

Due to the behaviour of L{X, a, x), it is quite clear that for any choice of .^ 7^ the func- 
tion 0,7 {X, a, x) presents a singular point, and therefore the potential Vc~{a,x) is always 
singular, in contradistinction to the case of V^{a,x) considered before. The presence 
of the singularity suggest that the results could be interpreted according to the method 
developed in p3[: the susy partner potentials are not directly related by isospectrality 
to the original potential V{a,x), but to a different problem consisting of this modified 
Poschl- Teller potential plus an infinite barrier potential placed precisely at the position 
where K~(a,x) has its singular point. 

12 



The case ^ = is interesting enough to be considered separately. It gives just the 
particular solution Vq~ {a, x), which is free of singularities. The susy partner potentials (pq) 



correspond to a factorization energy e~, which according to (25) is e~ = Eq. Therefore, 
there is an eigenstate ip oi H which is annihilated by the intertwining operator A^ ip{x) = 
0, and we can write 



4){x) tx exp 



Pq (y) dy 



(cosh ax 



il-A 



which is precisely the square integrable function ^0(2;) given in (^, the ground state of 
H. The eigenfunctions of Hq are then given by 

i^-{x) = iEn-e-)-'/^AQM^). n = l,2,..., (46) 

In the present case the possible missing state solving (|37|) would be 



if (x) (X exp 



Po (y) dy 



(cosh ax 



,A-1 



As A > 1, this function is not square integrable, and therefore it has not physical meaning 
as an eigenfunction of Hq with eigenvalue e~ . The spectrum of Hq is given simply by 
{En'-iU = 1,2, . . .}. Remark that, like in the previous case, this new Hamiltonian is not 
either strictly isospectral to H, although the reason is just the opposite: now the susy 
process eliminates one state of H without creating a new one which can substitute it, 
while in the previous situation a new state was created, but keeping the initial spectrum. 
Hence, we have another example of unbroken susy encoded in the spectrum of the couple 
H and Hq . 

It is interesting to evaluate the limit of the previous results for a ^ 00, in order to do 
that, we first analize the asymptotic behaviour of the factorization energy (|2J): 



-"-(Sf-i-^^- '-' 



For < a < g/A the susy potential Vq (a, x) in ( |45| ) is always attractive; on the other 
hand, when a > g/4: the potential becomes always repulsive; finally, when a = g/A 
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the potential vanishes identically. This change on the character atractive or repulsive of 
V(^{a,x) according to the values of a has a strong physical meaning. It is related to the 
fact that in the susy process, the ground state level Eq is eliminated from the spectrum 
of Hq , being always a member of the spectrum of H, irrespectively of the value of a. 
In particular, when a > g/A the potential V{a, x) has only this bound state (see coment 
after equation @), and therefore, in the same interval the new potential Vq {a, x) has not 
bound state at all. 

In Figure 3 we plotted three members of the original family (^ of modified Poschl- 
Teller potentials V[a,x) (the three thicker curves), and also, in the same type of lines, but 
thinner, their corresponding susy partner potentials V^~(a, x) given by (^). The values 
of the parameters are indicated in the caption. Remark that although the initial modified 
Poschl- Teller potentials are always negative, their susy potentials are less negative (in the 
case of the dotted line) or become even positive (as in the cases of dashed and solid curves). 

An important detail to be stressed is that in the limit a ^ oo the potential Vq (a, x) 
of ( p5[ ) has a well defined behaviour (unlike the situation for Vg (a,j;)): it becomes the 
delta barrier 

lim VQ{a,x) = +g5{x), g > 0. (48) 

It is quite remarkable the difference with the limit of the initial potential V{a,x), which 



was a delta well, as it has been shown in equation (11). 

Another interesting point to be considered is the analysis of the limit a ^ oo of (E^), 
which can be evaluated even for the case ^ 7^ 0. We already have all the information needed 
to write down this result, indeed: a{X — 1) ^ g/2, tanhax —>■ sgna;, plus equations (^), 
(H), and (H). We get the following: 

g e^l^l 
hm f3f (a,x) = - sgn x + f y — -- — r , (49) 

1- T7-r ^ ,, , , 2gge^Nsgnx , 2g^ e^^N 

hm V, (q, x)=g 6{x) + + -— — ^ -3 • (50) 



l-^(^)«gnx [l-C(^)sgn 
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Again, we have used the fact that the parameters a, A and g are related through @. Ob- 
serve that after the Umit process, we obtain a function with three different discontinuities: 

1. The function blows up like 2{x — Xs)~'^ at the singular point 



Xs = In 1 + — 

9 V Kl 

2. At the origin, we get a divergence of the type +g6{x). 

3. There is a finite jump discontinuity at the origin, due to the presence of the sign 
function in (50). This discontinuity is concealed by the presence of the superposed 



Dirac delta centred also at x = 0. 

These remaks are clearly illustrated on Figure 4, where we have plotted the limit case 
(without the delta distribution at the origin), plus one intermediate case. The dotted 
curve represents the plot of one of the susy partners of the modified Poschl- Teller potential 
(^5|), for the following values of the parameters: a = 1.9, A = 1.216 (or equivalently 
g = 1), ^ = —0.05. The solid curve represents the susy Dirac delta potential for g = 1 
and ^ = —0.05, and it is a limiting case of the dotted curve when a — > oo (the delta 
contribution comes out from the dotted hump). 

3.3 The connection with the susy Dirac delta potential 

Let us consider now the intertwining relationship (|l^) for the Hamiltonian associated with 
the Dirac delta well potential Vd{x) = —g5{x). Equations (|20| ) and ([2l| ) also hold, with 
Vd{x) as the known potential and the potential V to be determined. The relevant Riccati 
equation to be solved is now 

(3^-0 = -g6{x) - a, (51) 

where a is the factorization energy in this case. Remark that from the mathematical point 
of view it is a differential equation which includes a distribution. Therefore the solution 
could have some discontinuity. It is possible to find a particular solution of (pll), for a 
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particular value of the constant a, in terms of the sign function, indeed /3o(a;) = (5/2) sgn x 
for a = —((7/2)^. This function satisfies the differential equation almost everywhere, i.e., 
for every real value of x, except for x = 0. Then, the general solution can be found by 
using the standard technique of transforming the nonlinear Riccati differential equation 
into a linear one. The final result is the following: 

0oj(x) = — sgnx — In 1 — u; / e^'^' dy] = - sgnx H -, — —, — r . (52) 

^^ 2 ^ dx \ Jo V 2 ^ l_^(e£^^gg^^ ^ ^ 



But this is precisely the result previously obtained in ( |49| ) if we identify the two parameters 
uj = ^. Obviously, the potential coming out from this function will be exactly the same as 
(|50[). Therefore, we have been able to obtain also the general susy partners of the Dirac 
delta distribution as a byproduct of the general results derived for the modified Poschl- 
Teller potential in the previous section. We would like to insist on the fact that in the 
singular potential case (^ 7^ 0) the susy problem is related to a modification of the initial 
potential resulting from adding an infinite barrier placed at the singularity (see |23|). 

Let us comment now a little bit more on the results for the Dirac delta susy partners. 
On one side, the particular case obtained from ( ^0| ) by making ^^ = 

lim Vf^{a,x) = +g6{x) = V^'''^(x) (53) 

is in complete agreement with some already well known results y]. On the other side, 
the general solution (^ 7^ 0) presented in (|50| ) can be compared with a result recently 
published [^]. Some discrepancies can be appreciated between the results derived in this 
last paper and ours. Observe that the connection we have established between our results 
for the modified Poschl- Teller and Dirac delta potentials strongly supports the validity 
of our new susy partners for the delta. Finally, we would like to stress one of the main 
results: the susy partners and all the relevant information for the Dirac delta potential 
can be obtained taking the limit in the corresponding expressions for the Poschl- Teller 
case, due to the good behaviour of the limiting procedures. One remarkable difference is 
that for the delta potential there is just one susy partner (indeed, only one specific value 
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of the factorization energy a allows to obtain the solution of the Riccati equation (|5lD), in 
contrast to the analysis done for the modified Poschl- Teller potentials, where two different 
factorization energies e^ were found. 

4 2— Susy modified Poschl- Teller potential 

The higher order susy partners can be also determined for the potentials we considered 
before. Concerning this topic, it has been recently developed a handy technique using 
difference equations in order to construct multi-parametric families of isospectral potentials 



[ PQ , 21 1 . In this Section we shall comment briefly only on the results derived following 
this approach for the potentials we are dealing with in the 2-susy case. In the first susy 
step we use the factorization constant e^ and in the second step we use e~ (the same final 
second order susy results are obtained if the process is accomplished in the reverse order). 
Hence, the 2-susy potential is given by 

where V{a,x) is the initial modified Poschl- Teller potential (|I]), e are given by (^), and 
f3t{x), 137 {x) by (^) and (|4^), respectively. The 2-susy partner depends on the two 
parameters C, and ^. This family embraces a wide variety of potentials, one of the most 
interesting cases is obtained by taking only the particular solutions C = ^ = 0: 



ga 2a^ _ 2 ^(^ - 1) ^a^ 

2 cosh ax sinh ax cosh ax sinh ax 



Vb,o(a, X) = ;^ h ——2 = -a — -^ h ——^ • (55) 



This kind of solution is a particular case of a more general form of the modified Poschl- 
Teller potential used in some of the papers already mentioned |0] 



2/k(k + 1) A(A-1) 



y(^) = a' : ' - ' ' , (56) 

Vsmh ax cosh ax J 



precisely for the value k = 1. Note that (55) has not a well defined limit when a — > 00, and 
therefore it is not possible to find a 2-susy partner for the Dirac delta connected with a 
2-susy partner of the modified Poschl- Teller potential. Observe that this fact was implicit 
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from the begining, because only one factorization energy a was found for the 5, while for 
the modified Poschl- Teller we were able to find two different factorization energies e"*" , e~ . 

5 Final remarks 

Due to the relevance of solvable susy quantum mechanical models as toy examples for 
higher dimensional quantum field theories, and also for their use in solid state physics, we 
analized in detail the supersymmetry associated with the modified Poschl- Teller potential, 
which appears in many interesting physical situations, for example in the nonrelativistic 
limit of the sine-Gordon equation, in connection with a two-body force of Dirac delta type, 
when studying integrable many-body systems in one dimension, or when considering two- 
dimensional susy quantum field theories. 

In the present work we have constructed two new one-parametric families of exactly 
solvable potentials V^{a,x) and Vr{c(,x) related to the modified Poschl- Teller potential 
V{a,x) by the one-parametric superpotentials pf{a,x) and (37 {a, x), respectively. They 
represent two different cases of unbroken supersymmetry, and they reduce to some previ- 
ously published results for (^ = or .^ = 0. A relevant trait of our results is that, for specific 
values of C and a, the members of the family V/'{a, x) are free of singularities, which is, as 
far as we know, a fact unnoticed in the literature. On the other hand, the family Vr {a, x), 
with ^ ^ 0, embraces only singular potentials which have to be considered very carefully, 
because they are not susy partners of the initial potential (||) , but of the initial potential 



plus an infinite barrier at the singularity [23|. 



The connection between the modified Poschl- Teller and Dirac delta potentials was 
established, and we were able to construct the susy partner of the delta potential. The 
main remark is that only the singular family Vr{a,x) can be used to approximate an 
atractive delta potential in terms of the limiting procedure discussed in the paper. The 
explanation can be given in terms of the susy process: our particular solution of the 
Riccati equation ([Sl]) represents the superpotential f3o{x) = {g/2) sgnx usually derived for 
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the atractive delta potential y]. Therefore, the susy partner potential of the Dirac delta 
well is a Dirac delta barrier and, because this last potential has not bound states, the 
corresponding susy system has to present unbroken susy; in other words, the susy process 
eliminates the only bound state of the delta well in order to satisfy the Witten index 
condition for unbroken susy. The same holds for the case u; = .^ 7^ 0. The energy level at 
(T = Eq is destroyed and it does not play any role in the limit a ^> oo for the spectrum of 
VV~(a,x) (remember that linia^oo -E-q = Es). As the other energy levels dissapear, from 
the spectrum of V7~(q,x) when taking this limit, then the final potential has not bound 
states. The result is in complete agreement with the direct calculation. For V^{a,x), the 
situation is different because the ground energy level e'^ diverges as —a^ when a — > 00, 
and the potential has not a physically interesting limit. 
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2«r{A + l/2) 




Figure 1: The function M{X,a,x) given in Eq. pO|). Its 
bounded character enables the existence of non-singular susy 
partners V^{a,x). 



22 



-0.05 



-0.1 



-0.15 - 



E-i 




Figure 2: The susy partner potential V^{a, x) given by Eq. 
(^), with a = 0.1, A = 3, and ( = 0.0937, is an asymmetric 
double well. Its three bound states are also represented with 
dotted horizontal lines. 
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Figure 3: Different Poschl- Teller potentials V{a, x) from 
Eq. (H) (the three thicker curves), and their correspond- 



ing susy partner potentials Vf^{a,x) given by Eq. (45) (the 
three thiner curves). The values of the parameters are the 
following: dotted curves a = 1, A = 2.562, dashed curves 
a = 3, A = 1.758, and solid curves a = 6, A = 1.457. In all 
the cases the parameter ^ is taken to be zero. The initial 
Poschl- Teller potentials are always negative; their susy part- 
ners are less negative (dotted curve) or become even positive 
(dashed and solid curves). 
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Figure 4: The dotted curve shows the plot of a member of 
the family of susy partner Poschl- Teller potentials K~(a, x) 
given in Eq. (|45|), for the following values of the parameters: 
a = 1.9, A = 1.216 (or equivalently g = 1), and S, = —0.05. 
The solid curve represents the susy partner Dirac delta po- 
tential for g = 1 and ^ = —0.05, and it is obtained from the 
dotted curve when the limit a ^^ oo is considered. In the 
same plot there are two remarkable details: first the diver- 
gent term g5{x) (generated by the dotted hump) has not 
been represented, and second the potential has a disconti- 
nuity at X = 0, which is masked by the existence of the term 
g6{x). 
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